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Abstract: 

We present results on the two-loop leading and angular-dependent next-to-leading 
logarithmic virtual corrections to arbitrary processes at energies above the electroweak 
scale. In the 't Hooft-Feynman gauge the relevant Feynman diagrams involving soft and 
collinear gauge bosons 7, Z,W ± coupling to external legs are evaluated in the eikonal 
approximation in the region where all kinematical invariants are much larger than the 
electroweak scale. The logarithmic mass singularities are extracted from massive multi- 
scale loop integrals using the Sudakov method and alternatively the sector-decomposition 
method in the Feynman-parameter representation. The derivations are performed within 
the spontaneously broken phase of the electroweak theory, and the two-loop results are in 
agreement with the exponentiation prescriptions that have been proposed in the literature 
based on a symmetric SU(2) x U(l) theory matched with QED at the electroweak scale. 
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1 Introduction 



The main task of future colliders such as the LHC [E] or an e + e~ Linear Collider (LC) 
[ El El HI E| will be the investigation of the origin of electroweak symmetry breaking and 
the exploration of the limits of the Electroweak Standard Model. In order to disentangle 
effects of physics beyond the Standard Model, the inclusion of QCD and electroweak 
radiative corrections into the theoretical predictions is crucial. 

In the energy range of future colliders, i.e. at energies above the electroweak scale, 
y/s ^> M ~ Mw — Mz, the electroweak corrections are enhanced by large logarithmic 
contributions [ E] of the type 

These electroweak logarithmic corrections (EWLC) can be classified in a gauge-invariant 
way according to the powers iV of the logarithms of s/M 2 . The leading logarithms (LL), 
also known as Sudakov logarithms [Ej, correspond to N = 2L, the next-to- leading loga- 
rithms (NLL) to iV = 2L - 1, etc. 

The above logarithmic terms constitute the singular part of the corrections in the 
massless limit, M 2 /s — > 0. They result either as remnants of ultraviolet singularities 
after parameter renormalization, or as mass singularities from soft /collinear emission of 
virtual or real particles off initial or final-state particles. These latter do not cancel 
in observables, in contrast to the well-known soft and collinear singularities observed 
in QCD. This is, on the one hand, due to the fact that the masses of the weak gauge 
bosons provide a physical cut-off, and that there is no need to include real Z-boson and 
W-boson bremsstrahlung. On the other hand, the Bloch-Nordsieck theorem is violated 
also in inclusive quantities in non-abelian gauge theories if the asymptotic states carry 
non-abelian charges or in spontaneously broken abelian gauge theories if mass eigenstates 
result from mixing of gauge eigenstates [ E] • This leads to the appearance of LL also in 
inclusive quantities in such theories and thus in the electroweak Standard Model. As a 
consequence, the electroweak corrections can become of the order of the QCD corrections 
in the TeV energy range. 

These enhanced EWLC have found quite some interest recently; for reviews we refer to 
Refs. [E1E3- At the one-loop level the EWLC have been obtained, on the one hand, via 
explicit diagrammatic calculations for many 2 —>■ 2 scattering processes in the Standard 
Model and the Minimal Supersymmetric Standard Model [ EH E21 EH EH EH] • On the 
other hand, the universality of the electroweak LL and NLL has been proven, and general 
results have been given for arbitrary electroweak processes that are not mass-suppressed 
at high energies [EH1EZ! and applied, for instance, to gauge-boson pair production at the 
LHC [EH- 

The typical size of the one-loop electroweak corrections from LL and NLL for a 2 — > 2 
cross section is 

respectively, at yfs = 1 TeV, with M = Mw, and l — s 2 v = c 2 v = M^/M^. The size of the 
corrections increases with the number of particles in the final state, and it is important to 
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note that at the TeV scale, the LL and NLL have similar size and opposite sign resulting 
in large cancellations [ O EH] • 

Resummations of the EWLC have been proposed based on techniques and results 
known from QCD and QED. Fadin et al. [ 120] have resummed the LL by means of 
the infrared evolution equation (IREE), which describes the all-order leading-logarithmic 
dependence of a matrix element with respect to the transverse- momentum cut-off fi±, 
within a symmetric gauge theory This equation was applied to the electroweak theory 
by assuming that the /ij_-integration can be split into two regimes both corresponding to 
symmetric gauge theories. In the regime s > fi± > M, SU(2) x U(l) symmetry was used, 
whereas for M > fi± > A the U(l) em symmetry was assumed. Kiihn et al. have applied 
results from QCD to resum the logarithmic corrections to massless 4-fermion processes, 
e + e~ — > ff , up to the NLL [EI] and even to the next-to-next-to-leading logarithms (NNLL) 
[ I21j . This was done for a symmetric SU(2) x U(l) theory and additional electromagnetic 
effects were included following the IREE approach. It was found that at 1 TeV there is 
no clear hierarchy between LL, NLL, and NNLL. One of us has proposed a resummation 
of the NLL for arbitrary processes [ I22| 12*3*] . which relies on the prescription of matching 
a symmetric SU(2) x U(l) theory with QED at the electroweak scale. 

All these resummations amount to exponentiations of the EWLC. The approximate 
size of the two-loop LL and NLL resulting from the exponentiation of the one-loop cor- 
rections ()1.2|) for 2 — > 2 processes at y/s = 1 TeV is 

+ A 10 ** GsO = 3 - 5% - "S logS (sO = " 4 - 2% - (L3) 

respectively, and it is clear that in view of the precision objectives of a LC below the 
per-cent level these two-loop EWLC must be under control. 

All the above resummation prescriptions result from matching a symmetric SU(2) x 
U(l) theory and QED at the electroweak scale, and are based on the assumption that other 
effects related to spontaneous symmetry breaking may be neglected at high energies. This 
assumption needs to be checked by explicit diagrammatic two-loop calculations based on 
the electroweak Lagrangian, where all relevant effects related to spontaneous symmetry 
breaking are taken into account. In particular, the following non-trivial aspects need to 
be treated with care: (i) There is a multi-scale hierarchy of masses, M ^> frif^t ^ A, with 
heavy masses m t ~ Mh ~ M at the electroweak scale, light-fermion masses rrif^ t , and 
also an infinitesimal photon mass A, which is used as infrared regulator. As a consequence 
of this hierarchy, also logarithms of the large ratios M/mj^ t and mj^ t /A have to be taken 
into account, and the general form of logarithmic terms of order N in (jl.lj) becomes 
log^ 1 (s/M 2 ) \og N2 (M 2 /m 2 m ) log^ 3 (m^ t /A 2 ), with N = N x + N 2 + N 3 . (ii) In the gauge- 
boson sector, the gauge-group eigenstates B, W 3 mix resulting in mass eigenstates 7, Z 
with a large mass gap A M. (iii) Longitudinal gauge bosons appear as physical 
asymptotic states. 

The resummation of the LL has been checked for the massless fermionic singlet form 
factor in Refs. [ |2H I2S] and for arbitrary processes in the Coulomb gauge in Ref . [ |2H] • The 
resummation of the next-to-leading logarithms has so far not been confirmed by explicit 
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two-loop calculations. In the present paper, we consider a specific gauge-invariant subset 
of the next-to-leading EWLC to exclusive processes: the angular-dependent NLL of type 



where r represents a generic kinematic invariant different from s, and the ratio r/s depends 
on the angles between external momenta. These angular-dependent NLL are numerically 
important as has been stressed in Refs. [ElEHEj- Prescriptions for their resummation 
have been given in Refs. [EH OH] for massless 4-fermion processes and extended to arbi- 
trary processes in Ref. [ E3J - As we have mentioned above, these prescriptions are based 
on symmetric gauge theories. The purpose of this paper is to check them with an explicit 
two-loop calculation within the spontaneously broken electroweak theory. 

At one-loop order, in the 't Hooft-Feynman gauge, the angular-dependent NLL result 
only from diagrams where a gauge boson is exchanged between two external lines. Sim- 
ilarly, the angular-dependent NLL at two-loop order can be traced back to a relatively 
small set of Feynman diagrams. This allows us to present a diagrammatic calculation of 
the two-loop angular-dependent NLL for arbitrary processes. The calculation is based on 
the eikonal approximation. The relevant massive two-loop integrals are evaluated ana- 
lytically using two independent methods, one goes back to Sudakov [[7J, the other uses 
sector decomposition of Feynman-parameter integrals [ 123 I2H EH] • 

The paper is organized as follows: In Sect. EJ we define our conventions and the ap- 
proximations used in the high-energy limit and we discuss the Feynman diagrams that 
give rise to the leading mass singularities and the eikonal approximation. Section |3] is 
devoted to the description of the calculation of the two-loop integrals. The results for the 
contributing diagrams and their sum are presented in Sect. 0] The appendices provide 
information on our conventions and the explicit results of the individual loop integrals as 
well as relations between them. 

2 High-energy logarithmic approximation 
2.1 Preliminaries 

We consider generic electroweak processes involving n arbitrary polarized particles, 
which may be light or heavy chiral fermions, transverse or longitudinal gauge bosons, or 
Higgs bosons. As a convention, we consider n — > processes, 



where all particles (fi k and their momenta pu are assumed to be incoming. Corresponding 
2 — > n — 2 processes are easily obtained by crossing symmetry. Our calculations are 
performed in the physical basis, where the external particles (fi h as well as the virtual 
particles in the loops correspond to mass eigenstates, and mixing effects are properly 
taken into account. The matrix elements for the processes (J2.1j) and the external- leg 
gauge couplings are denoted with the shorthands 






(2.1) 
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where Ifq corresponds to the coupling of the gauge bosons V a = 7, Z, W 1 * 1 to an incoming 
particle ipi and an outgoing particle ipy. More details concerning gauge couplings can be 
found in App. |XJ 

All external- leg momenta are assumed to be on shell, p\ = mf , and we restrict ourselves 
to the high-energy region where s = {pi + P2) 2 ~ 2piP2 as well as all other kinematical 
invariants are much larger than the electroweak scale. In particular, we assume the 
following hierarchy of energy and mass scales 1 

\{p k + Pl ) 2 \ ~ \2p kPl \ > M 2 ~ Ml ~ > mj_, t > M 7 2 = A 2 . (2.3) 

The mass scale M is used to denote a generic weak-boson mass in the logarithms, and we 
neglect logarithms of the ratio M^/Mz, which originate from the difference between the 
Z- and the W-boson mass. With rrif^t we denote the masses of the light fermions. The 
infinitesimal photon mass A is used to regularize infrared divergences. 

All masses of real and virtual particles are assumed to be at or below the electroweak 
scale. Nevertheless, our results also apply to processes involving particles with masses 
that are heavier but of the same order as M, i.e. light Higgs bosons or top quarks with 
Mh — m t ~ M. However, the logarithms involving the ratios Mh/M and m t /M are 
neglected. 

Our next-to-leading logarithmic angular-dependent (NLL a ) approximation is defined 
as follows. We consider corrections that are logarithmically divergent in the limit where 
the ratios M/y/s, m/^t/M, and A/m^ t vanish. From these mass-singular logarithms we 
retain only the leading and the next-to-leading angular-dependent ones, i.e. at L loops 
(with L = 1,2) we consider only contributions of the order 



O 
O 



0<N<2L, and 

N / A/T2 



< N < 2L- 1, (2.4) 



where raii g ht,i are either light-fermion or photon masses. As stated above, the logarithms of 
ratios of heavy masses are neglected, i.e. we consider logMz — logMw — logMn — logm t . 

In our approximation all terms that are suppressed by factors M/s/s, m^/M or 
A/m/^t are neglected. In practice, all mass terms in the numerators of loop integrals are 
omitted in the calculations. In order to avoid factors of order y/s/M from longitudinal 
polarization vectors that would enhance mass-suppressed contributions, the Goldstone- 
Boson Equivalence Theorem (GBET) [ 30J has to be used for matrix elements involving 
longitudinal gauge bosons. For our purposes we can use the GBET in its naive lowest- 
order form since the quantum corrections to the GBET involve only two-point functions, 
which give no contribution in the considered NLL a order (|2.4j) in the 't Hooft-Feynman 
gauge. In practice each longitudinal gauge boson V{; = W L , Z L has to be substituted by 
a corresponding would-be Goldstone boson $ a = , \ using 

j^tpi 1 ...V£...<p in _ ^i-Q V a)j^<pi 1 ...^ a ...ip in ^ ^2.5) 



lr The first inequality implies, in particular, that all angles are larger than M/yfs in reference frames 
where all particle energies are not much larger than ^/i. 
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where Qv°- = ±1,0 is the corresponding charge. Thus, the general results for EWLC 
presented in the following have to be applied to the matrix elements involving would-be- 
Goldstone bosons. 

We restrict ourselves to matrix elements that are not mass-suppressed in the high- 
energy limit. This permits us to make use of the identity 2 



M 



Y^M' ■ I$ kik = O ( — J M' ll - in ~ 0, 



k=l 



(2.6) 



which can be derived from global SU(2) x U(l) symmetry and is very useful in order to 
simplify sums over external-leg insertions of the gauge-group generators. 



2.2 Feynman diagrams in eikonal approximation 

Mass singularities originate from diagrams with virtual particles coupling to on-shell 
external legs. In this paper we consider only NLL a contributions that result from the 
leading mass singularities, i.e. at two loops from contributions involving four mass-singular 
logarithms. We perform the calculation in the 't Hooft-Feynman gauge, where the gauge- 
boson propagators have the same pole structure as scalar propagators and the leading 
mass singularities originate only from diagrams with soft-collinear virtual gauge bosons 
coupling to different external particles. The exchange of soft-collinear scalar particles or 
fermions is mass-suppressed. The relevant two-loop diagrams have the structure 



U 2h,jk 



^ 3L, jkl 




D 



ab 

2C,jk 



r^abc 
u ?,Yjkl 




abc 
2Y, jk 



- Ly 4L,jfc«m 




where the soft-collinear gauge bosons V a , V b , V c = 7, Z, W 1 * 1 are exchanged between two, 
three, or four of the n on-shell external legs j,k,l,m = 1, . . . , n. The external lines, which 
are generically represented by full lines, can be fermions, transverse gauge bosons, would- 
be Goldstone bosons, or Higgs bosons. The external legs that do not couple to the soft 
gauge bosons are represented by the lines and the dots on the left-hand side of the grey 
blobs. 

In order to extract the leading mass singularities, the Feynman diagrams ()2.7|) are 
evaluated in eikonal approximation, i.e. by approximating the integrand in the limit where 
the momenta of the soft gauge bosons are small. In this limit, the above Feynman 
diagrams can be treated independently of the process and the spin of the external particles, 
i.e. universally for chiral fermions, transverse gauge bosons, or scalar particles. This is 
done as follows: 



2 Here the sums over the components i' k of the multiplets corresponding to the various external particles 
ik are implicitly understood. 
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• In the "hard part" of the diagrams corresponding to the grey blobs in ()2.7|) the 
momenta of the soft gauge bosons are neglected, and only the remaining "soft part" 
of the diagrams has to be integrated over the loop momenta. The blobs typically 
involve contributions from various tree diagrams, but they do not need to be eval- 
uated explicitly. In our derivations we only make use of the charge-conservation 
identity ()2.6j) to relate the complete tree-level amplitudes from different blobs. 

• The vertices involving three soft gauge bosons are associated with the usual Yang- 
Mills couplings 




= —iel"! 



C13 a 2 



9n 



-I 



(2.; 



where the particles and momenta are incoming, I^ a2 are the structure constants 
defined by (jA.3|) . and di indicates the complex conjugate of cij. In the 't Hooft- 
Feynman gauge, the propagators read — ig^ v /{I 2 — M 2 + is) for a gauge boson V a . 

• The emission of gauge bosons with momenta - along an incoming external 

line with momentum k\ = p ext in the soft limit l{ — *■ 0, 



fcl fc2 kz 

~* I ^ (2-9) 

gives rise to a product of terms containing a factor i/ (A;|— mf+ie) for each propagator 
with momentum kj = kj-i — and mass rrij and an eikonal factor [ IT7j 



Vi'i-kj+i) 



eik 



2k j iei^j, 



(2.10) 



for each vertex, where are the generators defined in App. EI Note that the 
form of these eikonal factors only depends on the gauge-group representation of the 
inflowing particles, but not on their spin. With = we denote equations that are 
valid within the eikonal approximation. 

• The eikonal factors defined above are proportional to the momenta of the emitting 
particles, and if these particles are virtual they depend on the loop momentum l\ 
via k 2 = Pe X t — h- If the integration is restricted to the region of soft gauge-boson 
momenta, l\ ~ 0, as in the Sudakov method, this loop-momentum dependence 
can be neglected, and one can use k 2 = p ex t in the eikonal factors. However, if 
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one integrates over the full loop-momentum space, as in the Feynman-parameter 
representation, also the region k 2 = p cy ± — l\ « where the emitting line becomes 
soft may be important. This happens for the diagrams D 2 q and -D3L, where such 
regions give rise to spurious leading logarithms if one uses eikonal factors with 
k 2 = Pext, whereas the loop- momentum-dependent eikonal factors defined in f!2.10|) 
suppress them 3 . The reason why these contributions have to be suppressed is the 
following: if the emitting line in (J2.9|) is a fermion or a scalar particle the mass 
singularity for k 2 = gets mass-suppressed in the complete Feynman diagram by 
the numerator. If the emitting line is a transverse gauge boson this remains true 
for the region k 2 ~ l 2 ~ 0, whereas the region with k 2 « and l 2 ~ p' eyLt 7^ p ext , 
which gives leading contributions only to the diagram D 2C , has to be suppressed in 
order to avoid double counting of topologically equivalent configurations when the 
sum over all soft-collinear gauge bosons is performed. 

• The denominators of the propagators denoted by full lines in (J2.7)) are kept exact, 
i.e. the square of the momenta of the soft gauge bosons is not neglected there. 

The explicit expressions for the diagrams ()2.7|) in eikonal approximation are given in 
App.|B21 

All terms involving four mass-singular logarithms originate from the diagrams (|2.7|) 
and there only from the terms that are kept in the eikonal approximation. Other con- 
tributions give rise to at most three mass-singular logarithms. In this paper we assume 
that all angular-dependent NLL result only from contributions with four mass-singular 
logarithms via the appearance of different scales in the logarithms. This is equivalent to 
the assumption that generic NLL are not multiplied by logarithms of ratios of kinematical 
variables that do not result from mass singularities. Although we have not proven this 
assumption so far, we do not see a source for additional angular-dependent NLL. We have 
checked for several examples that no such terms arise from neglected contributions to 
the diagrams (|2.7p . i.e. terms with gauge-boson momenta in the numerators. A complete 
proof of the assumption, however, requires a calculation of the complete set of NLL. 

2.3 Mass-gap effects 

Each topology in (|2.7|) has to be evaluated for all different mass assignments corre- 
sponding to the various electroweak gauge bosons V a , V b , V c = 7, Z, W*. In practice, for 
each diagram involving two or three soft-collinear gauge bosons we have the following 
four cases 

(M a , M b ) = (A, A), (A, M), (M, A), (M, M), (2.11) 
(M a ,M b ,M c ) = (M,M,M),(A,M,M),(M,A,M),(M,M,A), (2.12) 

whereas the external lines are assumed to have arbitrary masses m| at or below the 
electroweak scale. 

Our main aim is to investigate the effects related to symmetry breaking, and in par- 
ticular the effects of the large mass gap A < M in the gauge sector. This gives rise 

3 For all other topologies the loop-momentum dependence of the eikonal factors is irrelevant in NLL a 
approximation. 
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to logarithms of the photon mass A and light-fermion masses that violate SU(2) x U(l) 
symmetry, such that the full electroweak result is only symmetric with respect to the 
unbroken U(l) cm group. Nevertheless, according to the physical picture proposed in the 
IREE approach [ED] and generalized by other authors [HTH I2T1 (221 123]. the EWLC are 
expected to exhibit a higher degree of symmetry. In this picture, the complete electroweak 
result factorizes into a part which exhibits SU(2) x U(l) symmetry and corresponds to 
the full electroweak result for the case A = M, and a remaining part that originates from 
the mass gap A -C M and exhibits U(l) cm symmetry. 

In order to check this picture at the level of angular-dependent NLL, we organize our 
calculation as follows. All intermediate results /(A) depending on the photon mass A are 
split into two parts as 

/(A) = /(M)+A/(A), (2.13) 
where the part f(M) corresponds to the case A = M and has to be calculated for 4 

A = M > m m . (2.14) 

In this case, all mass singularities are regulated by M and the light fermion masses 
below the electroweak scale uif^t ^ M can be neglected. The remaining part, A/(A) = 
/(A) — /(M), originates from the mass gap A <C M. In the language of the IREE, this 
subtracted part can be understood as the part of the photonic contribution that originates 
from below the electroweak scale. 

2.4 Validity of our results for extensions of the electroweak theory 

Our derivations depend only on a few general features of the Electroweak Standard 
Model, such as the underlying global gauge symmetry, the spectrum of gauge bosons, 
and the fact that all particle masses are of the order of the electroweak scale or lighter. 
All leading and next-to-leading angular-dependent logarithms originate only from (soft- 
collinear) gauge bosons and depend only on gauge couplings. Therefore our results apply 
also to those extensions of the Electroweak Standard Model, where these features are 
preserved, i.e. where no additional gauge bosons and no logarithms involving mass scales 
much higher than the electroweak scale appear. 

Such models include softly broken supersymmetric extensions such as the Minimal Su- 
persymmetric Standard Model in the case where the masses of the superpartner particles 
are of the order of the electroweak scale M. Owing to mixing, the gauge couplings If fi 
may involve mixing matrices. More details about higher-order supersymmetric EWLC 
can be found in Ref. [I31j. 

3 Loop integrals in logarithmic approximation 

The loop integrals were evaluated in logarithmic approximation using two indepen- 
dent methods: the Sudakov technique and the sector-decomposition method described in 
Sect. 13 .11 and Sect. 13.21 respectively. These two methods were applied in a complementary 
way in order to calculate all integrals in App. [B]and to perform various cross checks. 

4 Note that this part f(M) cannot be obtained by simply substituting A = M in the result /(A) since 
the inequalities l|2.3(l and l|2.14|l exclude each other. 
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3.1 The Sudakov technique for angular-dependent logarithms 



The Sudakov technique has long been known and used for the calculation of the leading 
logarithmic asymptotics of field theories [ E2] . In this paper, we apply this approach 
to the calculation of the leading and next-to-leading angular-dependent logarithms at the 
two-loop level. To our knowledge, the Sudakov technique has so far only been applied to 
ladder or crossed ladder diagrams involving a single large invariant. Here we generalize 
this method to the case of different large invariants. 

We illustrate the method for the 3-leg ladder diagram 

Dstj ki show n in (ED) . In 

the eikonal approximation the corresponding integral is given by (|B.22|) . It involves four 
different mass singularities and thus gives rise to four large logarithms. These singularities 
appear if the gauge-boson momenta l\ and l 2 become soft and collinear to the momenta 
of the external particles to which the gauge bosons couple. Mass terms have to be only 
retained as far as they regularize the mass singularities. 

In order to extract these singularities it is convenient to use the following Sudakov 
parametrizations for the loop momenta, 

h = V\ \Pj - TT^-Pk) + x i [Pk ~ TT-^Pj] + K±, 
\ ZPjPk J \ ZPjPk J 

(TV?- \ ( tit? \ 

Pj ~ TT^Pl + x 2 [Pi ~ 7, ~Pj + h,±, (3-1) 
ZPjPi J \ tyjPi J 

where h,±Pj = h,±Pk — h,±Pj — h,±Pi = 0. The mass terms turn out to be only rele- 
vant for photon exchange since we assume rrij^i ~ M. The two-dimensional transverse 
momenta U t ± are space-like. They can be parametrized by their moduli |L ; j_| and az- 
imuthal angles 0j. Then up to irrelevant mass terms, the integration measures read 
d 4 /i = llpjpfcldl^d^idxid?/! and d% = l\pjPi\d^ _ L d02dx 2 dt/2- 

Since the leading logarithms originate from the regime of soft and collinear gauge- 
boson momenta, we can drop all /-dependent terms in the numerator and replace the 
gauge-boson propagators as 

> 7r5(2p j p k x 1 y 1 - M 2 a - l^ ± ), 
7r5(2p jPl x 2 y 2 -Ml-ll ± ), (3.2) 



ll - Ml + is 2p jPk x iyi -Ml-l 2 1± + is 



li - Ml + ie 2p jPl x 2 y 2 - Ml - ]* >± + ie 

up to irrelevant terms of order m^ kl . Performing the integrals over \f ± with the help of 
the 5 functions, we find after neglecting irrelevant mass terms 

4 r 
S 3h (M a ,M b ;pj,p k ,pi) = ^(PjPk)(PjPi)\PjPk\\PjPi\ J dxid|/idx 2 d?/2d0id02 

x 6(2p jPk x iyi - Ml) e(2p jVl x 2 y 2 - M 6 2 ) 

x [2 PjPkXl + m)y x + M 2 X\-^P 3 Pkyx - vn\x x + M^ 1 

x [2pjp k xx + 2pjpix 2 + rrijiyi + y 2 ) + Ml + Ml + 2p k pix x x 2 

+ 2p j p i x 2 y 1 + 2p j p k x 1 y 2 + 2l 1>± pix 2 + 21 2>± p k x 1 + 21 1> _ L l 2i± ] _1 
x [-2p jPl y 2 - m]x 2 + M b 2 ] -1 , (3.3) 
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where 1| ± = 2%\y\Pjp k — Ml and 1| ± = 2x 2 y2PjPi ~ M$ and the dependence on 0j 
enters via 1^. Terms involving four logarithms result from those parts of the integration 
region where all Sudakov variables are small, |xi|, |yi|, |x 2 |, \y 2 \ "C 1, and the integrand 
behaves as 1/ (£12/1X22/2), i-e. where each of the four denominators is dominated by a 
different term linear in one of the Sudakov variables. These parts of the integration 
region are selected by conditions of the form \2pjp k x\\ ^> |m|j/i|, \2pjpix 2 \ 3> \2pjp k xi\, 
etc. To leading-logarithmic accuracy, these conditions are implemented via step functions 
9(\2pjp k xi\ — |m|j/i|), 9(\2pjPix 2 \ — \2pjp k xi\), etc., which, in particular, ensure that none 
of the Sudakov variables can become zero. Upper integration limits are set to one, i.e. 
I X\ I, I yi I, I x 2 |, 1 2/2 1 < 1. Then irrelevant terms are neglected in the denominators and in 
the step functions. Since, in particular, all terms depending on and thus on 4>i are 
negligible, the integration over these angles can be performed trivially. After transforming 
regions with negative Sudakov variables to those with positive Sudakov variables one 
finally obtains (jB.23|) . 

In this approach it is crucial that the perpendicular components of the loop momenta 
\±_ can be dropped. This is the case for all ladder diagrams, if the Sudakov parametriza- 
tions are constructed from the external on-shell momenta of the lines to which the gauge 
bosons couple as in (|3.1j) . 

The Yang-Mills diagram S 2 y with two external lines contains three virtual gauge 
bosons that can become soft and collinear. If any two of them go on-shell, this results in 
four large logarithms. Thus, one has to sum over the three contributions with different 
pairs of on-shell gauge bosons to obtain the full leading- logarithmic result. Since there 
are only two relevant external momenta, the Sudakov parametrization is unique for each 
of the three cases. A slight complication arises from the fact that the numerator is linear 
in the loop momenta. Nevertheless, one can show that all terms involving transverse loop 
momenta can be neglected in the leading-logarithmic approximation and the leading- 
logarithmic contributions can be extracted in a way similar to the ladder diagrams. 

For the non-abelian graph S^y with three external lines no parametrization exists 
that would allow to neglect all 1$ j_ terms. For this diagram we therefore did not apply the 
Sudakov method but have checked the result from the sector-decomposition approach by 
a numerical integration of the Feynman-parameter integral. 

3.2 Sector-decomposition method 

The loop integrals were also evaluated in the Feynman-parameter representation. In 
this case, in order to extract the logarithmic mass singularities we used the sector de- 
composition [ 123 |2H1 l2Hj- which permits to factorize overlapping ultraviolet or mass 
singularities in Feynman-parameter integrals. A detailed description of this method is 
postponed to a forthcoming publication Here we only sketch the main steps of the 

sector-decomposition method applied to a generic two-loop massive integral with n + 1 
propagators 

r d% r d% Ngi&W) 

J {2-kYJ (27r)^n?i 1 (^ 2 -m? + k)' 1 ' ' 

where the momenta ki are linear combinations of the external momenta p\ and the loop 
momenta L, and the numerator iV is an arbitrary polynomial in these momenta. 
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Step 1: The integral is written in Feynman-parameter representation and split into 
n + 1 primary sectors as described in Ref. [I2H]- After eliminating the usual S function 



5(1 -E 



n+1 



Xi j 



each primary sector gives rise to a Feynman-parameter integral over the 



unit cube in n dimensions of the form 



[o,i]" [D{x)X 



D(x) 



sP s (x) + rPJx) + ... + m 2 P m (x) + \ 2 P x (x) 



(3.5) 



where the resulting denominator D(x) is split into polynomials according to the hierarchy 
of scales in the diagram, which is assumed to be s > r > . . . > m 2 ^> A 2 in (J3.5J) . Note 
that in order to extract the angular- dependent logarithms log (s/r) with r = t, u, . . ., we 
compute the integrals in the Euclidean region in various limits of the type s ^> t = u, 
s = t 3> u, etc., where we separate the energy scales in various ways. 

Step 2: The polynomials in (j3.5J) have various zeros if subsets of Feynman parameters 
vanish, e.g. 



P(x) = 0, at x\ 



x a 



0. 



(3.6) 



which give rise to mass singularities. In the presence of such a singularity, the polynomial 
can be written as 



P{x) = Y J x i P i {x). 



(3.7) 



In order to separate the singularity associated to (j3.6j) from other overlapping singularities, 
we decompose the corresponding integration domain [0, l] q into q subsectors Qj with 



> x^j, and in each subsector Vtj we perform variable transformations x^ 



Xjx[ for 



all i 7^ j, which remap £lj — > [0, l] q and bring the polynomial (J3.7|) into the form 



P(x) 



where the variable Xj is factorized. 



x 



x 



(3i 



Step 3: Recursive application of step 2 permits to factorize all zeros at all scales, 
until the denominator assumes the form 5 



D(x) 



sP s (x) Y[ x/ + rP r (x) 



+ + rn 2 P m (x) A x? + X 2 P x (x) 



k=l 



1=1 



(3.9) 



where 0^,6^,0^ are positive integers. In (|3.9|) all Feynman parameters that give rise to 
mass singularities are factorized and the polynomials P are non-vanishing. This allows 
for a simple power counting of the logarithmically divergent integrations. 

Step 4: All logarithms of ratios of scales can now be extracted in NLL a approximation 
by analytical integration, where the polynomials P can be treated as constants P(x) ~ 
P(0). At present, explicit results are available [EH] for the special class of integrals where 

5 In general the denominator assumes the form D(x) = D(x) Yii=i x i , an d the overall factorized 
Feynman parameters Xi with di > can be cancelled by corresponding terms in the numerator. 
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the various subsets of parameters {xj\aj > 0}, {xk\bk > 0}, etc. that are associated 
to different mass scales are disjoint. As an example, for the case of a hierarchy of four 
different scales s ^> r ^> M 2 ^> A 2 , we have integrals of the type 6 



d l x / d m y / d n z 



s e 






e-1 


{ 


(s U l i= i Xi + r) 




n^i^ + A 2 } 6 



n n+m—p i i 
NLL a ^2 



1 



p=0 5=0 



p\ q\ (N-p-q)\ l0gP I A 2 / "'^ V .} [- J 



(3.10) 



where iV = /+m+ra and e > 1. Such integrals permitted us to calculate all diagrams listed 
in App. IB. 21 except for the ladder diagrams with simultaneous photon- mass and external- 
mass singularities, since these diagrams lead to integrals where the various subsets of 
parameters {xj\aj > 0}, {xk\bk > 0}, etc. are not disjoint. Such integrals have not been 
solved so far, since it was more convenient to perform the calculation using the Sudakov 
method. 



4 Results 



In this section we present results for one- and two-loop Feynman diagrams evaluated 
in the high-energy limit ()2.3|) using the eikonal approximation (eik) and to next-to-leading 
logarithmic angular-dependent (NLL a ) accuracy ()2.4j) . All results are split according to 
()2.13|) into contributions corresponding to A = M and remaining A contributions, which 
originate from the gap A C M in the gauge sector. 

In Sect. I4.ll we first recall the one- loop results [US], we then present in Sect. 14.21 our 
results for various subsets of two-loop diagrams and for the complete two-loop corrections, 
and in Sect. 14.31 we discuss the exponentiation of these logarithmic corrections. Explicit 
results for the individual one- and two-loop integrals can be found in App. El 



4.1 One-loop results 

Within the 't Hooft-Feynman gauge, the one-loop LL and angular-dependent NLL [ 
ITT)] originate from diagrams where a gauge boson V a = 7, Z,W ± is exchanged between 
two different on-shell external legs j ^ k, 




In the eikonal approximation, these yield 

SM im = ^.M;/ (4.2) 

a 

6 Here we consider all terms with the total power of logarithms equal to N = I + m + n, but in NLL a 
approximation we only need the contributions with N — p — q < 1. 
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where the integral S(M a ;pj,p k ) defined in (|B.2|) can be decomposed into 



S(M a ; Pj ,p k ) N ^ La E(M a ; rrij) + E{M a - m k ) + R(M a ; Pj ,p k ), 



(4.3) 



with 



A 



A 



M 2 



E(X; rrij) = - log 2 log 2 E(M; rrij) = E(M) = - log 2 



in 



M 2 s 
R(M a ; Pj ,p k ) = 2 log — £ log 



(4.4) 



The functions E depend only on the energy scale s, on the internal masses M a , and on 
the masses rrij ik of the external lines, and do not give rise to correlations between the 
two external legs j, k in the sum (14.5)1 . whereas the function R contains logarithms of 
PjPk/s depending on the angle between the momenta pj and p k , but is independent of the 
external masses. 

The complete one-loop correction is obtained by taking the sum over all pairs of 
external legs, and the part originating from the functions E can be simplified by using 
the charge-conservation identity ()2.6|) . This yields 



5M1 1 - 



1 n 

2 £ SM \m 



a 

47T 



-i=i 



i 1 ...i....i„ 



t> x; Mir 

3,k=l 



C^E(M) + 5 i?j Ql AE(X; rrij) 



£ Ih.Il ih R(M; Pj ,p k ) 

a=^,Z,W 3 



+ ^,'; 0,,^[ ;,Qi, AR(X;pj,p k ) 



(4.5) 



where C ew represents the electroweak Casimir operator defined in (|A.5|) . Qi is the eigen 



value of the charge operator H. 



and 



AE(X; = E(X; mj) - E(M; mj) = log 



'mf 



log 



M 2 



X 2 \ 



AR(X;pj,p k ) = R(X;pj,p k ) - R(M;pj,p k ) = 2 log y— I log 



\ 2 PjPk\ 



(4.6) 



In order to discuss the two-loop results in the next section, it is useful to rewrite the 
one-loop result ()4.5|1 in matrix form. To this end we introduce the following notation 



M = M 11 - 



MI a (k) =M h - 



.i k ...i n ja 

etc., 



MI a (k)I b (k) = M* 1 " 4 *'"^!^!^, 



(4.7) 



where the generators I a (k) are matrices acting on the indices corresponding to the external 
leg k of the matrix element, with commutation relations 

~I a (k),I b (l)]=5 kl £ P{k)P cb . (4.8) 
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Using this notation also for other matrices as C ew (k) or Q(k), the result ()4.5|) can be 
rewritten as 



SMi = Mo Sew = M (<5 sew + 



(4.9) 



where the complete electroweak (EW) result is split into a symmetric electroweak (sew) 
part 



"sew 



ewIa=m 



(4.10) 



£ { - J E ^ ew 0') ^s 2 — + E E ^(jW) io t 



M 



47T 



i,fc=i a=7,Z,W± 



3 — 



which corresponds to the case A = M and is manifestly SU(2) x U(l) symmetric, and a 
remaining subtracted electromagnetic (sem) part 



A(5 EW = Sew - Sew\ x =m = T~ { ~ \ E ^0) 



47T 



2 f=i 



m 2 A 2 , 2 m 2 n 
2l0g V l0g M^" l0g M 2 



+ E Q(j')Q(^)iog 



s A^ 
|2p iPfc | ° g M2 



(4.11) 



which originates from the gap A <C M, i.e. from the A terms in ()4.5|) . and is U(l) em 
symmetric. 

4.2 Two-loop results 

In the following we present results for the two-loop diagrams (12. 7j) combined into three 
subsets where the soft-collinear gauge bosons couple to two, three, or four of the n on-shell 
external lines, respectively. 

The two-loop results are decomposed into reducible contributions, which involve prod- 
ucts of the one-loop integrals (|4.3jl . plus remaining irreducible parts. This decomposition 
is based on the explicit two-loop results in the NLL a approximation ()2.4|) given in App. IB. 21 
and the relations given in App. O 



Terms from two external lines 

We begin by considering the contributions 



a,b 

corresponding to the diagrams 
)b (a + 



SM'ftft = E D&jk + D&jk + E [D&jk + D 



abc 
2Y,kj 



a,b 




k 




(4.12) 




(4.13) 
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i.e. diagrams with n on-shell external legs and soft-collinear gauge bosons V a , V , V c 
7, Z, W exchanged between only two of these external lines. These yield 



5M m) = U 1 M ° 



(\ \ 2 . .i\...i'....i' k ...i r , 



a,b \ 3 3 k k 



b]Pj,Pk) 



+ E 



(7*1=),, ijL IZ c S 2Y (M a ,M b ,M c ;p 3 ,p k ) + (j <- k) 

j J k 



V47T/ 



a, 6 ^ K 



a 3 



i ™ 

+ -V V T a 

* a,c h=i 

hjtj,k 



I^Jl lk I^AS 2C (M, X; Pj ,p k ) + (j <- fc) 



(4.14) 



where we have made use of the identities ^M, (EH3j) . (JOJ), (jHljl . (RT2|) . and (|Tl3|) . 

Terms from three external lines 

Here we consider the contributions 



tMisft = E E ^« + E r > 

a,6 ^7r(j,fc,0 c 



(4.15) 



where we sum over all six permutations 7t(j, k, I) of external lines j, k, I. These contribu- 
tions correspond to the diagrams 



E 

u,b 



k 



+ (j <-> k) + (j <->/)+£ 




(4.16) 



with exchange of soft-collinear gauge bosons V a , V b , V c = 7, Z, W between three external 
on-shell lines, which yield 



x vf" ••" — ( —\ 2 \a : ■■■■■>', ->'i. ■■■>',-'•■■ 



E E ( /6ja L. J |i fc 4^3L(M a , M 6 ; 



Pj,pk,pl) 



. 7r(i,fe,0 a,6 J J 

+ Y,Ih ] IbJi l Jl H SMM a ,M b ,M c ;p J ,p k ,p l ] 
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An J 



7rO',A,0 ^ a,b i 3 

+ E ( /a/7 ),',, IhJlS{M-p v p k )^S{X-p^ Pl ) 



- \ E AS 3L (M, ^'iPhPk, Pi) } , (4-17) 



2 

^ a,c 

where we used f|P3|L (fOjl . (Q, 0, O, and (fC~7jl 

Terms from four external lines 

Finally, we have the contributions 




sM&iL) = E D&jum = E IT f ' ( 4 - 18 



l 2(jr'fcim) 

a, 6 

originating from gauge bosons coupling to four external legs, which reduce according to 
(IB.32|) to simple products of one-loop integrals 

a, b 

(4.19) 



Complete two-loop correction 

We now combine the contributions from the above subsets of diagrams into the com- 
plete virtual two-loop correction to an arbitrary process involving n on-shell external legs 
as follows 

in i n i n 

5M n-.in = Y!M%$ + - E'-^feS + g E' M%% m) , (4.20) 

j,k j,k,l j,k,l,m 

where we have to sum over all combinations of external legs with appropriate symmetry 
factors. The primes indicate that the sums include only terms with different external legs, 
i.e. 

n n n n n n 

E' == E , E' == E , E' == E ■ (4-21) 

j,k j,*=i j,k,l j,k,l=l j,k,l,rn j,k,l,m=l 

k^j kj^j;l^j,k k^j;l^j,k;mjtl,j,k 

We first consider the irreducible contributions to ([4.20)1 . i.e. the contributions from 
A^c in ()4.14|) and ASsl in (|4.17|) . which could not be expressed as products of one- loop 
integrals. These contributions cancel, 
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x [AS 2C (M,X; Pj ,p k ) - AS 3L (M, X; Pj , P i, Pk )} N = a 



(4.22) 



because of (|C8|) . The complete two-loop correction is thus given by the reducible contri- 
butions to ()4.14j) . ()4.17|) . and ()4.19|) . i.e. contributions from products of one-loop integrals 
S. These yield 



5M\ 1 - in N = a 



x S(M a ;p j ,p k )S(M b ;pi,p r , 



Era ja jb jb 



j,k,l 



a.b 



3 3 



+ E ( ia n 9j %^(M; pj,p h )AS(\;pj,pi) 



+ EM."' 



a.b 



3 J 



+ ( /S/7 V, 7 (^ 7 W S(M; Pi> p fc ) AS(X; Pj ,p k ) 



NLL, 



1 / a 

2 Vl^J 



i , > «i...«:...ij c ...« n 



EX 



E ( /a/S )i'i, E ( M ) + ( J7/7 )i'i, A£(A;m,; 



x 



£(/»/*).,. £(M) + (rr%, ifc A£(A;m fc ) 



6 ' 



+ 2 £ (rrPP)^ E(M)AE(X\ m,) + (PPPP),,. [A£(A; m^f 

a 

l/a\ 2 " 

) X>o J " ' . 

b 



2 V4tt 



E(/^L ^W + (^ 7 )^ A£(A;m,) 
h i 1 



x 



E q^RiMwM) + IlJl k AR(X; Pj ,p k ) 



3 l k lk 



(\ \ 2 , J i l ...i'j...i' k ...i n 



An 



em 



E( ja/6j5 kv% ^( M ) W 



a, 6 



+ Y,( lbli n^% k E{M) AR(X; Pj ,p k ) 

a 

+ E( /a/7/7 k^ fc A ^( A ; ™;) R(M; Pj , Pk ) 



+ (PPP)^!^ AE(X; rm) AR(X- Pj ,p k ) 



(4.23) 
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where we have used f!2.6|) . ()4.3|) . ()A.3|) . ()A.4|) . and the fact that R(M a ;pj,p k ) is symmetric 
with respect to interchanging j «-> fc. Note that the photonic coupling matrices J 7 / 7 are 
to the right of I a in the term containing AE(X; rrij) R(M\pj,pk). 

Using 7J. = —Qi^i'i- and rewriting ()4.23|) in the compact operator form introduced 

j 3 i 3 

in (|4.7|) . results in 



£M 2 ^ _.M (^ cw + 25 scw 4 cm + <&mj , (4.24) 

where <5 sew and 5 sem are the operators defined in (j4.10|) and ()4.11|) and correspond to the 
symmetric electroweak and the subtracted electromagnetic parts of the one-loop correc- 
tions. It is important to note that these two operators have a non-vanishing commutator 



|2P*P||, l-N s ». M 



iog^^iog 3 - w — niog 



rn 



2 

light, i . 



N = l,2, (4.25) 



where mug^t is either a light-fermion or a photon mass. This means that the ordering 
of the terms 5 sew S sem in ()4.24|) . which is determined by the contribution involving I a PP 
in (|4.23p . is relevant at the level of angular- dependent NLL. As we point out in the 
next section, the determination of this ordering constitutes an important aspect of our 
result and permits to discriminate between different exponentiation prescriptions for the 
electroweak corrections. 



4.3 Exponentiation 

If we combine the two-loop correction (|4.24|) with the one-loop correction ()4.9|) and 
the Born amplitude, to two-loop NLL a accuracy we find the exponentiated form 

M 2 = Mo + 5Mt + 5M 2 N = a M exp (5 SCW ) exp (5 scm ) . (4.26) 

In particular, the form of the two- loop correction operator (S^ ew + 2S sew S sem + S^ em ) /2 
implies that the symmetric electroweak part 5 sew and the subtracted electromagnetic part 
5 sem exponentiate separately, and that the latter exponential is external. This means that 
the charge operators in exp(5 sem ) can be identified with the charge eigenvalues of the 
external particles in the process. 

At the level of the LL, this result confirms the exponentiation of the EWLC obtained 
with the IREE [ 2(1 and already checked for arbitrary processes by a two- loop calculation 
in the Coulomb gauge [ I26j . We found also agreement with the results of Refs. [ 12*4*1 123] 
for the special case of the fermionic form factor corresponding to the decay g — > ff of 
an SU(2) x U(l) singlet g into massless fermions. We have explicitly verified all results 
of Refs. [|2U 123 by evaluating the subset of diagrams _D 2 l, D 2 c, and D 2 y in (|2.7|) for the 
special case of massless external particles. At the level of angular-dependent NLL, our 
result is in agreement with the exponentiation prescriptions adopted in Refs. [ El |^] for 
massless fermionic processes, and extended in Ref. to arbitrary processes. 

This agreement indicates that, at least up to the level of angular-dependent NLL, 
a symmetric SU(2) x U(l) gauge theory matched with QED at the electroweak scale 
provides a correct physical picture for the resummation of EWLC in the high-energy 
limit. This picture has been formulated within the theoretical framework of the IREE 
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[ !2Uj. which describes the all-order leading-logarithmic dependence of matrix elements 
on the transverse- momentum cut-off fi±. This infrared scale fi± is the crucial ingredient 
in order to avoid the difficulties related to the breaking of the SU(2) x U(l) symmetry 
that originate from the large mass gap A < M in the gauge-boson sector. In fact, the 
scale fi± permits to separate two regimes of the electroweak theory both with exact gauge 
symmetry. The regime \/s > fi± > M, which is insensitive to the gauge-boson masses 
and has SU(2) x U(l) symmetry, and the regime M > /i^ > A, where the weak gauge 
bosons are "frozen out" and only U(l) em symmetry is left. 

To our knowledge, the IREE has been formulated only at the level of LL, and the 
application of the physical picture described above to the level of NLL relies on a weaker 
theoretical basis. At this level, the following two arguments can be used for the exponen- 
tiation of the next-to-leading logarithmic corrections. On the one hand, if A = M then the 
SU(2) x U(l) symmetry is restored in the gauge sector and one expects the exponentiation 



as in a symmetric SU(2) x U(l) theory. This permits to predict the two-loop term pro- 
portional to S 2 ew in ()4.26|) and implies that 5 scm = at A = M. On the other hand, 
the logarithms of the photon mass and light-fermion masses originate only from photons 
coupling to external legs and are expected to exponentiate as in QED. However, the QED 
results can be generalized to the electroweak corrections only if the contributions from vir- 
tual photons can be separated from those of the weak gauge bosons in a gauge-invariant 
way. This is the case only if s = M 2 , where the logarithms of s/M 2 originating from 
virtual weak bosons vanish. Here one expects 



where 5 sem corresponds to the QED corrections. This, together with (j4.27J) . permits 
to determine <5 sem and the two- loop term proportional to <5 2 em in ()4.26|) . However, we 
note that the above two conditions, (j4.27J) and (j4.28J) . are not sufficient to determine 
the two-loop interference terms S sew 5 sem between symmetric electroweak and subtracted 
electromagnetic contributions, which vanish in both cases A = M and s = M 2 . These 
two-loop interference terms are an important result of our electroweak calculation for 
s 3> M 2 ^> A 2 . In particular, they are crucial in order to predict the ordering of the two 
exponentials in ()4.26|) . which starts to be non-trivial at the level of angular-dependent 
NLL as indicated by the commutator ()4.25|) . 

5 Conclusions 

We have studied the two-loop asymptotic behaviour of virtual electroweak corrections 
to arbitrary processes involving light or heavy chiral fermions, transverse or longitudinal 
gauge bosons, or Higgs bosons. We have calculated the two-loop leading and angular- 
dependent next-to-leading logarithmic contributions in a process-independent way in the 
region where all kinematic invariants are much larger than the electroweak scale. The 
relevant Feynman diagrams involving exchanges of soft and collinear virtual gauge bosons 
7, Z, and W between on-shell external legs have been evaluated in the eikonal approxi- 
mation in the 't Hooft-Feynman gauge. 



A = M M 2 = M exp [S sew ] , 



(4.27) 



s = M 2 ^ M 2 = M exp [<5 sem ] 



(4.28) 
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Analytical expressions for the relevant two-loop integrals, which involve up two six 
different scales, have been calculated via two independent methods. On the one hand, 
we have evaluated the integrals in the Feynman-parameter representation using sector 
decomposition to isolate the mass singularities in the integrand and performing the inte- 
gration in logarithmic approximation. This method was applied to all diagrams except 
for those ladder diagrams with simultaneous photon-mass and external-mass singulari- 
ties. On the other hand, we have employed the well-known Sudakov method, which is 
very efficient for calculating diagrams with only one large energy scale but turns out 
to be more complicated for diagrams with more large energy scales. In particular, in 
the Sudakov approximation we did not succeed in calculating the diagram where three 
soft gauge bosons interacting via a Yang-Mills vertex couple to three different external 
legs. In all diagrams where both methods could be applied we found agreement at the 
angular-dependent next-to-leading logarithmic level. 

In order to isolate the effects originating from the large mass-gap between the photon 
mass A and the weak-boson masses M w ~ M z ~ M, which breaks the symmetry in 
the gauge-boson sector, the loop contributions depending on the photon mass have been 
split into a part corresponding to A = M and a remaining subtracted part. Combining 
the results from all diagrams we found that the sum of the two-loop leading and angular- 
dependent next-to-leading logarithmic corrections can be written as the second-order term 
of a product of two exponential functions. The first exponential contains the part of the 
corrections corresponding to A = M, i.e. the SU(2) x U(l) symmetric part. The second, 
outer exponential contains the contributions that originate from the mass gap A <C M 
and corresponds to the QED corrections subtracted in such a way that they vanish at 
A = M. 

This result agrees with resummation prescriptions that have been proposed in the 
literature. These prescriptions are based on the assumption that, in the high-energy limit, 
the electroweak theory can be described by a symmetric SU(2) x U(l) theory matched 
with QED at the electroweak scale, and that no additional effects from spontaneous 
symmetry breaking appear. Our result, which has been derived within the spontaneously 
broken phase of the electroweak theory and in the physical basis, demonstrates that this 
assumption is correct at the next-to-leading angular-dependent logarithmic level. 

Our derivations depend only on a few general features of the Electroweak Standard 
Model, i.e. on the underlying gauge symmetry and the fact that all particle masses are 
of the order of the electroweak scale or lighter. Therefore our result is also valid for 
those extensions of the Electroweak Standard Model that contain only novel particles 
with masses of the order of the weak scale and no additional gauge bosons. Such models 
include, for instance, the Electroweak Standard Model with two Higgs doublets or softly 
broken supersymmetric extensions such as the Minimal Supersymmetric Standard Model 
in the case where the masses of the Higgs bosons and the superpartner particles are of 
the order of the electroweak scale. 
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A Gauge-group generators 

All our derivations are performed in terms of the physical (mass-eigenstate) gauge 
bosons 7, Z,W ± . The corresponding gauge couplings result from combinations of the 
generators T a and Y of the electroweak SU(2) x U(l) gauge group, and read 

r = _q, f * = i«-4q t^ = ±^ = ±^. t (a.i) 

where Q = T 3 + Y/2 represents the electric charge, and we use the shorthands c w = cos# w 
and s w = sin# w for the weak mixing angle, which is fixed by = 1 — = M^/M| in 
the on-shell renormalization scheme. 

For the matrix components of the generators (jA.l|) we use the notation 7^, where 
a = 7, Z, W ± denotes the gauge fields, whereas the indices %' and % correspond to two 
physical (mass-eigenstate) components </?£/ and ipi of a multiplet. The explicit matrix 
representations corresponding to the scalar doublet, right- or left-handed fermions and 
gauge bosons, as well as more details concerning our conventions, can be found in App. B 
of Ref. [H71. 

The matrix component 1^ determines the gauge coupling for the vertex with the 
particles V a and <£>j incoming and the particle </v outgoing. The fields and the matrix 
components are in general complex and satisfy the relations [ IT7] 

fe)* 'I',- ~I% (A.2) 

where the particles a, i, . . . correspond to the charge conjugated of a,i, . . .. 
In our derivations we make extensive use of the commutation relations 



ja rt 



E /C/ c fe (A.3) 

c=7,Z,VK± 



and of the well-known commutation relations 



£ [l a f a , I b ] = with b = 7, Z, W ± , (A.4) 

=y,Z,W ± 



for the electroweak Casimir operator 



1 (Y 



C™:= ^ iafa = -rU) +^ T (T+l\ (A.5) 

where T is the total isospin, and T(T + 1) is the Casimir operator of the SU(2) group. 
The electroweak Casimir operator is a diagonal matrix apart from the neutral gauge-boson 
sector, where mixing gives rise to the non-diagonal components CJZ = Cf^ = —2c w /s w 
(see App. B of Ref. [Ej). 
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B 1- and 2-loop integrals in logarithmic approximation 

In this section we present detailed results for the one- and two-loop integrals involving 
soft-collinear gauge bosons. For each diagram we first specify the corresponding Feynman 
integral in eikonal approximation (eik.). Then we also give the corresponding integral in 
the Sudakov approximation (Sud.). Finally, we present explicit results in next-to-leading 
logarithmic angular-dependent (NLL a ) approximation ([2.4)1 . which have been obtained in 
the high-energy limit ()2.3|) and for all cases specified in (j2.11j) and (j2.12|) . These results 
were derived using the Sudakov approximation and the sector-decomposition method 
described in Sect. I3.ll and Sect. 13.21 respectively. 

The external momenta are assumed to be on-shell, pi = ml, with masses at or below 
the electroweak scale, i.e. M > m k ^> A. Masses that do not regularize mass singularities 
are neglected. Consequently, the masses of the external particles and of the internal 
particles that are not soft-collinear gauge bosons are only relevant for photon exchange 
diagrams, where the masses before and after photon emission are equal. Therefore, we 
can set the internal and external masses of the particle lines equal in the following. 

B.l One-loop integrals 

For the one-loop diagram 




(B.l) 



we have the integral 

Here and in the following we suppress the infinitesimal imaginary parts ie of the causal 
propagators for brevity. In the Sudakov approximation 

S(M a , Pj ,p k ) S ± d (B.3) 
_ f 1 dxi r 1 dyi ( M 2 a \( m) \ / m\ \ 

Jo x x Jo yi \ iyi \2pjp k \) V 12^*1 7 \ 1 \2pjPkl 1 )' 

The results in NLL a approximation corresponding to the cases M a = M, A, are given in 
(JSJ) and (tOt . 
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B.2 Two-loop integrals 

2-leg ladder diagram Szl 

We begin with the planar ladder diagram 

■3 



n ab 




V47T/ 

x Mt " 



M b ;pj,p k ) 



with 



d 4 /i fd% 



S 2h (M a ,M b] p v p k ) := -(AttYJ—-LJ- 



(2tt)V (2^ [^-M^fe-ZO 2 -^] 



[(p fc + hf - mm ~ Mi][(Pj -h- h) 2 - m|][(p fc + h + hf - m^] 
In the Sudakov approximation 

S 2L (M a ,M b ;pj,p k ) u " 

= 4 



1 dxi /- 1 di/i /- 1 dx 2 Z" 1 dy 2 



o xi Jo ?/i Jo x 2 Jo y 2 



M, 



\2PjPk\J 



x 2 y 2 



\2PjPk\) 



M 2 



x 0(x 2 - x l )6(y 2 -y l )d[x 1 - 



m) 



\1PjPk\ 



yi f 2/i 



x hr 2 



'/7J 



^JfeM te- 



rm 



x 2 i 



In NLL a approximation, we find the following expressions for the cases (j2.11j) : 



NLU 1 i 2 A 



5 2L (A, X; Pj , Pk ) ™ - log 2 |2 ^ fc| _ ( ^ )2 6 _ |2w| 



x 



21og 3 ^ J * +3 log 



log" 



m 2 m 2 1 



log 



A 2 



rrijTnl 



5 



log — — —r + log 



(IPjPk) 2 
2 



^ 77T. 2 

log 2 — — — + log 



\2PjPk\ 



\ 2 PjPk\) 



2 ™fc 



ni 



\ 2 PjPk\ 

2 

log- 



mi 



+ \og z 

,2 



2 



log 2 



mi 



6 |2pj-p fc | |2pj-p fc | 



i 2 m 7 

log — — ^— + log 



S 2L (\,M; Pj ,p k ) 



NLL a 



6 l0g |2pj-p fc | 3 



|2PiPfe| 
- ir log 3 



\2PjPk\ \2PjPk\ 

2 V 



2 ™fc 



|2pjPfe|/ 



1 , , M 2 2 , , M 2 



log 



2 2 

m^mi 



\2pjp k \ (2pjp k y 



2 ° g |2 W | 



r?j 



log r— 3 i +log 

|2piPfc| 



mt 



|2pjPfc| 



-2 log 



A 2 



log 



2 2 

m1jm k 



\2pjp k \ (2pjp k y 



NLLa a »/ \ NLL a 1 1 4 ^ 



5 2L (M,A; Pj ,p fe ) = a S 2h (M,M; Pj ,p k ) = a - log 4 ^ 
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2-leg crossed ladder diagram S 2 c 

For the 2-leg crossed (non-planar) ladder diagram 



D 



2C,jk 



we have 

S 2C (M a ,M b ;pj,p k ) 

x 




) 2 S 2C (M a , 



M b ;pj,pk) 



x Mir'' ■■'""(I''r )r,,(r !'■),■„ 



(B.10) 



d% r dH 



{2<kYJ {2-KY[ll-Ml][{ Pj -l l f-m 
l6\Pj(Pk ~ h)]\pk(Pj ~ h)] 



[H ~ M?] [(p* - l 2 y - ml] [fo + l 2 f - mf ] [{jp k + h - l 2 ) 2 - mfj ' 



;b.h) 



As discussed in Sect. I2.2[ in order to avoid spurious leading logarithms originating from 
the region l\ pa pj and l 2 pa p k when the integral is evaluated in the Feynman-parameter 
representation, loop-momentum-dependent eikonal factors (I2.10J) have to be used for the 
inner vertices of this topology. 

In the Sudakov approximation, where the loop-momentum dependence of the eikonal 
factors can be neglected, we have 



S 2 c(M a ,M b ;p j ,p k ) = 

ridxi r^-dyt f 1 dx 2 r 1 dy 2 ( M 2 \ ( 

= 4/ / / / S iiji-r: 1 \9\x 2 y 2 - 

Jo xi Jo ?/i Jo x 2 Jo y 2 \ \2pjPk\J \ 



Mi \ 



\ 2 PjPk\J 



x9[yi 



2 



xi\6\x 2 



ml 



rV2 0(y 2 -yi)9(x 1 -x 2 ). 



\2pjp k \ ) \ \2pjPk\' 
In NLL a approximation, we derive the following expressions for the cases (|2.11|) : 

a i \ \ \ NLL a 1 , ^ 2 , 3 m ] m \ 

5^2c(A, \;Pj } pk) = - log , „ = i log 



;B.12) 



3 \2pjp k \ (2pjPk¥ 



+ log 4 



24 V 6 \2p jPk \ 

5 lo ^ lo 

6 ° g \2pjp k \ ° g \2pjp k \ 



\ 2 PjPk\J 



\ 3 



- - log 



2 m j i 2 

log 



mi 



\ 2 PjPk\ 



4 \2pjp k \ \2pjp k \ 
+ log 2 



log 2 ^ ■ ml 



S 2C (\,M;pj,p k ) = S 2C (M, X;pk,Pj, 



NLL a 



log : 



\ 2 PjPk\ 

M 2 



log- 



m 2 



3 |2pjp fc | |2pjp fc | 



NLLa 1 i 4 M 



(B.13) 
(B.14) 
(B.15) 
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2-leg Yang-Mills diagram S 2Y 

For the 2-leg Yang-Mills diagram 



D 



abc 

2Y, jk ■- 




cik / a 



v — j S 2Y (M a ,M b ,M c ; Pj ,p k ) (B.16) 



we have 



S 2Y (M a ,M b ,M c ; Pj ,p k ) := (4vr) 4 1 



d 4 /x /■ d 4 / 



(2^7 (2 7 r)4[/f-M2][(p,-/ 1 ) 2 -m J 2 
'PjPk[(k + 2h)Pj] 



[11 - Mi][(h + l 2 y - M?][(Pj + h? - mMPk - k) 2 - ml 



(B.17) 



For the cases ()2.12|) we find the Sudakov approximation 



Sud 



S 2Y (M a ,M b ,M c ; Pj ,p k ) 

rl dx\ r 1 dyi r 1 dx 2 f 1 dy 2 



Jo xi Jo y 1 Jo x 2 Jo y 2 \ 



VTLa 



\1pjpk\ 



ml 



X6 \ V2 1 2pjp k 1 2pjp k | 



V2 



x 9 x 2 y 2 



M. 



\ 2 PjPk\J 
Ml 



x 2 y 2 



Ml ' 
\2PjPkl, 



0(Vi ~ 2/2) +d[x 1 y 1 - 



ill 



Ml 



\2vM) 6{x '~ Xl)e{m ~ Vl) ~^ xm ~Wk) 



\ 2 PjPk\J 



Ml 



x 9 x 2 y 2 - r 9 x x y 2 - 



Ml 



Ml 



Ml 



\ 2 PjPk\J 

and in NLL a approximation we obtain 

S 2Y (X,M,M; Pj ,p k ) N ^ La S 2Y (M,X,M;p k , Pj ) 



\2pjp k \ \2pjp k \ \2pjp k \ 



NLL a 



S 2Y (M,M,X;pj,p k ) S 2Y (M,M,M; Pj ,p k ) 



1. , M 2 



log J 



log 



|2pjPfc| |2pjp fc |' 

NLL, 1 



■ ]oe 4 M 2 
6 |2pjp fc | 



(B.18) 



(B.19) 
(B.20) 
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3-leg ladder diagram 

For the 3-leg ladder diagram 



D 



3L, jkl 




S 3L (M a ,M b ;p j ,p k ,pi) 

...l' r ..l' k ...i' v ..ln a ( T bja\ jb 



(B.21) 



we have 



S 3h (M a ,M b ; Pj ,p k , Pl ) := -(47r) 4 | J 



16(PjPk)[(Pj + h)pi] 



{( Pk - hY - mm - M 2 ][( Pl + l 2 y - m]][( Pl - hY - mf] ' 



B.22) 



As discussed in Sect. 12.21 in order to avoid spurious leading logarithms originating from 
the region l\ ph —pj and l 2 ~ when the integral is evaluated in the Feynman-parameter 
representation, a loop-momentum-dependent eikonal factor (|2.1(Jj) has to be used for the 
emission of the gauge boson V b along the line j in this topology. 

In the Sudakov approximation, where the loop-momentum dependence of the eikonal 
factor can be neglected, we have 



S^M^M^p^p^pi) = 

Jo X\ Jo yi Jo x 2 Jo U2 \ 



M. 



n 



\2PjPk\J 



x 2 y 2 



Ml \ 
\lPiPi\) 



xe (x a -^4x l )e(x 1 



m A 



x9[yi 



\PjPi\ 



mi 



\2PjPk\ 



xi)0[y 2 



\2PjPk\ 
\ 2 PjPi\ 



Vi )0 \x 2 



mi, 



\2PjPi 



■2/2 



x 2 



;B.23) 



Neglecting angular-dependent NNLL of order log 2 (s/M 2 ) log 2 (2p m p n /s), we obtain 
the following results for the cases ()2.11|) in NLL a approximation: 

c ( \ x % NLL a 1 i m ) m l , m^mf A 2 

^3l(A, X;Pj,Pk,pi) = ^log^^log^log ^ 



1 m 2 mf, A 2 

~2 og (2^F og WM 



i 2 m l 
log — — —r + log 



2 m f 



\ 2 PjPi\ 



\ 2 PjPi\ 



+ log 



2 2 

m-mi 



log- 



\PjPi\J. 



1 

+ r 



i 2 m i i 2 m l 

log 2 J + log 2 



l 2 PiPi| 



2 

log 2 — — — - + log 2 



m 2 k 



+ log 



\^PjPk\ L3 



log 



mi 



\^PjPk\ 
\PjPk\~ 



\2PjPk\ 



mi 



— 1 

12 ° g |^| J 
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+ log : 



+ log 



2 m ) 



\2pjPk\ 
\2p&k\ 



3 , 2 m k 3 , o m i 
log 2 fc - - log 2 



9 i„„ i_ m > 



log 



log 



8 |2pjp fe | 8 \2pjp k \ 4 Ipjp/I |2pj-p & | 



1 i 3 m k 
* log 3 



1 , 3 m / 
-log 3 



6 \2pjPk\ 6 |2p^| 



mr 



log 



lo \PjPk\ / I 1q 
1 



+ -lo 



+ log 2 



log 



mr 



|2PiPfc|) 



m, 1 , 2 m i i 2 m 
log 2 log 2 



24 |2pjp fe | 6 |2pjp fc | |2pjp fe | 8 \2pjp k \ \2pjPk\ 



+ log 



\pjPk\ (1,3 ™i 



log 2 



n "<7 x i 3 

log — -i-r - - log 3 



2 \ 

i \ 



+ (m fc (^) 2 - MPiPk?) W ^ + ^ log g log 3 ^ j , (B.24) 



S 3L (M,X;pj,p k ,pi) 



NLL, 



log 



+ 77I0 



M 2 



3 |2pj-pi| 6 |2pjPi|J |2p 3 -pi| 



log 3 



M 2 



(B.25) 



S 3L (X,M;p j ,p k ,pi) N = a log 1 log 



2 2 
rrijiJik 



log 2 



M 2 



1 4 M 2 
7 l°g 



1 2 M 2 
-2 log 



l 2 PjPi| 



[2pjP*| (2pjPk) 2 \2pjPi\ 6 |2pjp fc | 



2 m ? 

log — — r + log 



m 



M 2 



|2pjPfc| 



\ 2 PjPk\J 



log' 4 log 1 — - 

3 |2pj-pjfc| |2pjPfc| 



S 3h (M,M;pj,p k ,pi) 



NLL a 1, .4 M 2 



2 bg ISpiftl 



(B.26) 
(B.27) 



3-leg Yang-Mills diagram S 3Y 

For the 3-leg Yang-Mills diagram 



r^abc 




we have 



S 3Y (M a ,M b ,M c ; Pj ,pk,pi) := (An) 4 J 



—J S 3Y (M a ,M b ,M c ; Pj ,pk,pi) (B.28) 

...i^...iy..i n j}y jq j a 

X ^0 1 i t .i j 1 i' k i k 1 i' l ii 1 bc 



d 4 /l /■ d 4 /; 



(2tt) 4 J (27r) 4 [/ 2 -M 2 ][(p J -/ 1 ) 2 -m 2 ] 
8pfePi[(ii + 2/ 2 )Pj] - 8pjPk[(h - - 8p i p i [(2/i + / 2 )Pfe] 



R 2 2 - Mi][(h + l 2 f - M?][(Pk - h) 2 - m 2 k ][(pi + h + k) 2 - m 
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(B.29) 



Neglecting angular- dependent NNLL of order log 2 (s/M 2 ) log 2 (2p m p n /s), we obtain the 
following results for the cases (|2.12|) in NLL a approximation: 



5 3 y(A, M, M;pj,p k ,pi) = S 3Y (M, \,M;pi,pj,p k ) = S 3Y (M, M, X;Pk,Pi,Pj) 
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NLL a 1 \PjPk\ , 2 M 
= -o lo g^^ lo S To 1 

3 \pjpi\ \2p k pi\ 



M 2 m 2 
log — 1 — 3 log : 



\2PkPi\ 



\2PkPl\ 



S 3Y (M, M,M] Pj , Pk , Pl ] 



NLL a 



0. 



(B.30) 



4-leg ladder diagram S 4 l 

Finally, for the 4-leg ladder diagram 




k eik / a 



x M 



S 4 h{M a , M b ;pj,p k ,pi,p 



m ) i 



(B.31) 



"j—"k—"l—"m — '"n. ja ja 



ln ja ja jb jb 



we have 



Pj > Pk j Pi j Pm 

) = S(M a ; Pj ,p k )S(M b ;pi, Pm ). 



(B.32) 



C Relations between loop integrals in logarithmic approximation 

In the following we list the relations between the one- and two-loop integrals that have 
been used in Sect. 14.21 in order to simplify the sum over all eikonal contributions. These 
relations have been obtained from the results of App. IB. II and App. IB. 21 and are valid in 
NLL a approximation 7 ()2.4j) . 

Combinations of 2-leg ladder integrals can be expressed as products of one-loop inte- 
grals using 



[SzL(M a , M b ;p j7 p k ) + S 2C (M a , M b ;pj,p k )] + (a <-> b) N = a S(M a ;p j ,p k )S(M b -p j ,p k ), 

(C.l) 

which is valid for all cases (j2.11j) . The 2-leg Yang-Mills diagram can be related to the 
2-leg crossed ladder diagram using 



S 2Y (M aj M bl M c] p v p k ) N = a -^S 2C (M a ,M b ; Pj ,p k ) 



(C.2) 



which is valid in all cases (12.121). Furthermore, we have 



S 2L (M,X; Pj ,p k ) N = a S 2L (M,M; Pj ,p k ) 
S 2C (M,\;pj,p k ) = S 2 c(\M;p k ,p j ). 



(C.3) 



The relation 



5 3L (M a ,M 6 ;p i ,p fe ,p z ) + S 3L (M b ,M a ;pj,p h p k ) N = a S(M a ; p j: p k )S(M b ] p^ p,) (C.4) 



7 The relations are actually valid for the 0-function representations given in App. [B] as well, which 
contain also NNLL. 
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permits to simplify combinations of 3- leg ladder integrals for all cases ([2.110 . Furthermore 



£ sgn(7r(j,k,l))S 3h (M,M; Pn ,p^, P7r .) N ^ 0, (C.5) 

w(j,k,l) 

where the sum runs over all permutations ir(j, k, I) of j, k, I, and sgn(7r(j, k, /)) is the sign 
of the permutation. For the 3-leg Yang-Mills diagram we have 

S 3Y (M, M, M; p 3 ,p k , Pl ) N ^= La 0, (C.6) 

as is evident from the totally antisymmetric property in the external momenta Pj,Pk and 
P\. In presence of a photon, this feature is absent owing to the mass gap A <C M, and we 
find instead 



Ssy(M, A, M;pj,pk,pi) N = a - 



A5 , 3L (M, \;pj,pi,Pk) - AS 3L (M, \;pi,Pj,Pk) 



(C.7) 



where AS3L is the subtracted part of the 3-leg ladder diagram as defined in (j2.13[) . This 
is related to the subtracted part of the 2-leg crossed ladder diagram by 



ASsl (M, A; Pj,Pu Vk) N = a AS 2C (M,\;p vPk ). (C. 
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